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We revisit the statistical significance of the "dark flow" presented in Kashlinsky et al. (20091. We 
do not find a statistically significant detection of a bulk flow. Instead we tind that CMB correlations 
between the 8 WMAP channels used in this analysis decrease the inferred significance of the detection 
to 0.7a. 

Subject headings: cosmology: observations 



1. INTRODUCTION 
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A recent set of papers (Kashlinsky et al. 
claims to have detected the velocities of galaxy clusters 
with respect to the cosmic microwave background (CMB) 
frame by means of the kinetic Suny aev-Zel'dovich (kSZ) 
effect (Sunyaev & Zeldovich 1972). The papers sug- 
gest the existence of a "dark flow": a 700 km s _1 bulk 
flow of all matter out to a redshift of at least z ~ 0.1 
(r ~ 400 Mpc) . The magnitude and direction of the flow 
are claimed to be consistent with the peculiar velocity 
of the Local Group with respe ct to the CMB fram e as 
inferred from the CMB dipole (Kogut et al. 1993[ ). Ve- 
locity coherence over such large scales is not predicted by 
the standard ACDM cosmology and would, if confirmed, 
constitute a major observational result. 

In this paper we revisit the analysis presented in Kash- 
linsky et al.f ( |2009[ ), hereafter referred to as K09. The K09 
analysis seeks to measure the kSZ signal of a sample of 
~700 X-ray-selected galaxy clusters. The 3-yea r WMAP 
temperatu re maps for 8 differencing assemblies (Hinshaw 
et al. 2007 1 are high-pass filtered in an attempt to remove 
the primary CMB anisotropy. The temperatures of the 
filtered maps at the galaxy cluster locations are fit to a 
dipole, which is interpreted as the kSZ signal induced by 
a bulk flow of the galaxy clusters. 

We will argue that the uncertainty of this measurement 
is dominated by primary CMB anisotropy, not detector 
noise. As the CMB is observed by all 8 WMAP channels, 
the errors are highly correlated between these channels, 
and the inferred detection significance is greatly reduced. 

2. CLUSTER SAMPLE 

We construct a cluster sample as similar to that used 
in K09 as possible. For all clusters we require that 
z < 0.3 and a corrected X-ray flux in the 0.1-2.4 keV 
band > 3 * 10~ 1 ^ergsJ 1 cm" 2 . We use the REFLEX cat- 



alog (Bohringer et al. 2004) and require that 5 < 2.5°, 
leaving 415 clust ers. We "use the eBCS catalog ( |Ebeling 
et al.||1998j |2000| ) and require that 5 > 0° and \ b\ > 20", 
leaving 279 clusters. We use the CIZA catalog ( |Ebeling 
Kocevski et al.|[2007 l, which, after our cuts, 



et al. 2002 



contains 122 clusters. Mote that this versio n of the CIZA 
catalo g contains 13 clusters total (73 from|Ebeling et al. 
( |2002"t and 57 from |Kocevski et al.| ( |2007| )), whereas K09 
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used an extended version containing 165 clusters total 
which is not publicly available at the time of this writing. 
We note that while K09 removed all clusters whose X-ray 
emission appeared to be dominated by a point source, we 
do not . This does not have a strong effect on our best-fit 
dipole, which, as discussed in Section [4] is very close to 
the value presented in K09. We use 816 clusters total, 
compared to 782 in K09. 720 of our clusters survive the 
3-year WMAP KP0 galactic mask used in this analysis. 

3. CMB MAPS 

Our preparation of the filtered CMB maps duplicates 
the method used in K09. We use the "foreground re- 
duced" t emperature maps from the 3-year WMAP data 
release 1 ( |Hinshaw et al.||2007| . We use one map each 
from the 8 differencing assembly channels: Ql, Q2, VI, 
V2, Wl, W2, W3, and W4. For each channel we con- 
struct the filter Ft described in K09. The filter is es- 
sentially a high-pass filter with a transition multipolc of 
I ~ 300. We explicitly remove the I = 0, 1, 2, 3 compo- 
nents from each map, as described in Kashlinsky et al.] 
(2008), by setting F(l= 0,1, 2, 3) = 0. We use the 3-year 
WMAP KP0 galactic mask for all steps of the analysis 
which involve spherical harm onic transforms. W e use the 
HEALPix software package (Gorski et al. 2005). 

The next step is to construct a mask which isolates the 
temperature fluctuations at the cluster locations. This 
mask has all pixels outside of the cluster areas and out- 
side of the KP0 galactic mask set to zero. The "clus- 
ter areas" are defined in K09 to be disc-shaped regions 
surrounding each cluster with rdisk =min[60x-ray, 30'], 
where 6*x- ra y is related to each cluster's X-ray emission. 
Because we lack access to the list of #x-ray used in K09, 
we use rdisk = 30' for all clusters. As pointed out in K09, 
because the majority of clusters have ^x-ray > 5', the 
K09 analysis uses rdisk = 30' for most clusters. Specifi- 
cally, the average rdisk used in K09 is 28'. 4, the standard 
deviation is 3'. 2, and only 16 clusters have rdisk < 20'. As 
such, our method is very similar to the method used in 
K09. As discussed in Section|4] any difference in method- 
ology has little impact on our best-fit dipole, which is 
very similar to that presented in K09. 

The final step is to calculate the dipole component 
of the temperature fluctuations at the cluster locations 

1 http: / /lambda. gsfc.nasa.gov 



2 



(which have been isolated using the pixel mask de- 
scribed in the previous paragraph). This is accom- 
plished with the HEALPix function remove-dipole, which 
returns a 3-component dipole vector [a x ,a y ,a z ]. The 
dipolc is calculated for each of the 8 WMAP chan- 
nels individually and results from all channels are com- 
bined using inverse-variance weighting. For example, 

a s = (E,i w x,iax,i)/(J2i w x,i)> where w x>i = cr^j and i 
is the index for the 8 WMAP channels. The variances 
are calculated from simulations which are described in 
Section |5] 

4. RESULTS 

Our best-fit dipole is [a x ,a y ,a z ]=[1.2, -2.4, 0.2] pK. 
This may be compared to [a x ,a y ,a z ]=[0.G, -2.7, 0.6] /iK, 
the corresponding all-z result from K09. The magnitude 
(direction) of our best-fit dipole is within 6% (17°) of 
the best-fit dipole measured in K09. We conclude that 
our best-fit dipole agrees well with that measured in K09, 
suggesting that the slight differences in methodology (the 
cluster catalogs and the choice of r<ii B k) are unimportant. 
Furthermore, we have repeated this analysis with rdisk 
increased and decreased by 20%, and the results do not 
change significantly. The magnitude (direction) of the 
best-fit dipole changes by less than 4% (11°) compared 
to the r^k = 30' case. 

Our best-fit dipole may be described as a vector with 
2.7/iK magnitude and with a higher temperature in the 
direction of the galactic coordinates (£,b) — (298,4). A 
naive interpretation of this dipole suggests a bulk flow 
moving away from (£, b) — (298, 4), which is the opposite 
sign of the velocity presented in K09. However, the in- 
terpretation of the sign and magnitude of the velocity is 
complicated by the convolution of the kernel of Fg with 
the gas profiles assumed for the clusters. This issue was 
not discussed in K09. For these reasons we do not at- 
tempt to quantify the sign or magnitude of the inferred 
velocity. We feel that this is justified by the low detection 
significance of the dipole, as discussed in Section [7] 

5. ERROR ESTIMATION 

We estimate the error of the dipole measurement as 
follows. The basic strategy is to repeat the analysis de- 
scribed in Section [3] but with simulated WMAP data 
replacing actual WMAP data. We generate 1000 real- 
izations of the CMB using the best-fit 3-year WMAP 
ACDM Cf T power spectrum. We use 3-year results, as 
opposed to 5-year, to remain consistent with the anal- 
ysis of K09. We convolve each CMB realization with 
the beams of the 8 WMAP channels in order to sim- 
ulate noise- free observations. For each CMB realization 
we also generate white noise maps for each WMAP chan- 
nel. The noise maps are generated using the prescription 
outlin ed in the WMAP T hree- Year Explanatory Supple- 
ment (Limon et al. 20081, with the noise variance in a 
given pixel inversely proportional to the number of times 
that pixel was observed. Our simulated noise is un- 
corrected between WMAP channels and between map 
pixels. While the latter assumption is not strictly true 
for WMAP data, the effects are negli gible for any tem- 
perature analysis ( Limon et al. 2008 1 . The white noise 
maps are added to the simulated noise-free CMB maps 
to produce maps which should have the same statistical 



properties as the actual "foreground reduced" temper- 
ature maps described in Section [3] We have confirmed 
that the (KPO-masked) power spectra of these simulated 
maps match those of the true maps. 

These simulated maps are then analyzed using the 
methods described in Section [3] To summarize, the maps 
are filtered by Fg and the dipole component of the filtered 
maps at the locations of the galaxy clusters is calculated. 
Although the WMAP data is simulated, we use the ac- 
tual cluster positions and the actual KP0 galactic mask. 
We note that because the definition of Fg given in K09 
depends on the measured spectrum C| fey , we calculate 
Fg for each simulated map. 

This suite of 1000 simulated WMAP "experiments" 
allows us to measure the uncertainty with which each 
WMAP channel measures each dipole component. As 
there is only "noise" (CMB and detector noise) in these 
simulations, the distribution of dipole measurements pro- 
vides a measure of the uncertainty of a single measure- 
ment of the dipole. These distributions are well de- 
scribed as Gaussian with zero mean. Inverse-variance 
weighting is used to combine the different channels' es- 
timates of each dipole component. Finally we are left 
with 1000 measurements of the 3 dipole components. 
These distributions are also well described as Gaussian 
with zero mean and are shown in Figure [l] The la 
measurement error for each dipole component is defined 
to be the best-fit Gaussian width a of each distribu- 
tion. The final uncertainties on the dipole components 
are [c 0a .,cr a , <7 a J=[1.7, 1.7, 1.1] pK. These uncertainties 
should be accurate to within a few percent. The uncer- 
tainty is highest on the dipole components that lie in the 
galactic plane (a x and a y ) because of the geometry of the 
galactic mask. 

6. CMB CORRELATIONS 

These simulations bring to light an important fact that 
was not discussed in K09: the dipole estimates are highly 
correlated between the 8 different WMAP channels. The 
correlation coefficient p between two different channels' 
estimates of a given dipole component is approximately 
0.9. For example, the correlation coefficient p between 
the estimates of a x provided by the Ql and W2 channels 
is 0.90, as shown in Figure [2j We attribute these correla- 
tions to primary CMB fluctuations, as the detector noise 
is uncorrelated between channels. 

To test the hypothesis that these correlations are 
caused by CMB fluctuations, we have separately fil- 
tered the noise-free CMB maps and the detector-noise- 
only maps in these 1000 experiments. The filter Ft is 
still constructed using the full maps which include the 
CMB and detector noise. The resulting uncertainties 
on the dipole components in the CMB-only maps are 
[a ai , a ay , cr a J=[1.7, 1.7, 1.1] pK. These uncertainties are 
consistent with those obtained using the full simulations 
which contain both CMB and detector noise. The uncer- 
tainties on the dipole components in the noise-only maps 
are much smaller: \a am , a ay , cr a J=[0.2, 0.2, 0.1] pK. We 
conclude that the uncertainty of the dipole measurement 
is dominated by CMB fluctuations, not detector noise. 
Although the filter Ft is constructed with the intent of fil- 
tering out primary CMB anisotropy, residual CMB power 
is still present in the filtered maps and dominates the er- 
ror on the dipole measurement. 
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Furthermore, if CMB fluctuations dominate the un- 
certainty of the dipole measurement, then the 5-year 
WMAP maps should produce a best-fit dipole that is 
very similar to that obtained using the 3-year WMAP 
maps. We have repeated ou r analysis on 5-year WMAP 
maps 1 ( |Hinshaw et al.| |2009) to test this hypothesis. We 
use the 3-year KPU mask and 3-year Fi filters on the 
5-year maps in order to make the comparison as direct 
as possible. The best-fit dipole is [a x ,a y ,a z ]—[1.3, -2.3, 
0.1] /xK, which is very close to the best-fit dipole from 
the 3-year WMAP data. The magnitude (direction) of 
this dipole is within 2% (4°) of that obtained using the 
3-year WMAP maps. This provides further support to 
the claim that the uncertainty of the dipole measurement 
is dominated by CMB fluctuations, not detector noise. 

If we simulate 1000 WMAP "experiments" and en- 
force the unphysical condition that each WMAP chan- 
nel observes a different realization of the CMB, then the 
dipole estimates are uncorrelated, as expected. In this 
unrealistic scenario the uncertainty on the dipole mea- 
surement is much smaller: [a aw , cr a , er a J=[0.7, 0.7, 0.4] 
/xK. These errors are closer to those presented in K09: 
[0a x , <7 ay , c a J = [0.5, 0.4, 0.4] /xK. These errors are smaller 
than the errors described above (which account for CMB 
correlations) by a factor of ~ v / 8 ; as is expected when 
combining 8 uncorrelated estimates as opposed to com- 
bining 8 highly correlated estimates. 

7. DETECTION SIGNIFICANCE 
Our best-fit dipole is [a x ,a y ,a z ] — [1.2 ± 1.7, —2.4 ± 
1.7, 0.2 ±1.1] fiK. The x 2 / d -0-f- is 2.52/3. The probabil- 
ity to exceed this x 2 is 0.47, corresponding to a Gaussian 
detection significance of 0.7c. If we use the best-fit dipole 
presented in K09, [a x ,a y ,a z ] = [0.Q, -2.7, 0.6] /xK, which 
is quite close to ours, the detection significance is 0.8cr. 
The errors used in these significance calculations come 
from the simulations presented in Section [5] which take 
into account correlations between the WMAP channels. 



A slightly different statistic may be considered if one 
is specifically interested in constraining bulk flows: the 
component of the best-fit dipole projected along the di- 
rection of the peculiar velocity of the Local Group with 
respect to the CMB frame . Measurements of the CMB 
dipole (Kogut et al. 1993) suggest that this velocity is 
towards the galactic coordinates (i,b) — (276,30). This 
statistic has 1 degree of freedom and we have calculated 
the uncertainty on its measurement using the methods 
described in Section [5] The best-fit projected dipole is 
2.2 ±1.6 /xK, corresponding to a detection significance of 
1.4ct. The sign is such that the temperature is higher at 
(£,b) = (276,30). 

We conclude that there is not a significant detection 
of a bulk flow. The significance of the best-fit bulk flow 
is 0.7a and the significance of the component projected 
along the Local Group's peculiar velocity is 1.4o\ 

8. CONCLUSION 



We have revisite d the analysis presented in Kashlinsky 
et al.| ( |2008[ |2009[ ) which reports a significant detection 
of a bulk flow ot ~700 galaxy clusters out to z ~ 0.1 
by means of the kSZ effect. We have demonstrated that 
the estimates for the kSZ signal are highly correlated 
between the different WMAP channels used in this anal- 
ysis and that this correlation is caused by primary CMB 
anisotropy. We have simulated the errors on the kSZ 
measurement while taking into account these CMB cor- 
relations and find that there is not a significant detection 
of a kSZ signal or bulk flow. 
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Fig. 1. — 1000 simulated estimates for the 3 dipole components. These simulations take into account the CMB correlations between the 
different WMAP channels. The uncertainty is highest on the dipole components that lie in the galactic plane (a x and a y ) because of the 
geometry of the galactic mask. 




Fig. 2. — 1000 simulated estimates for the a x dipole component from two randomly chosen WMAP channels, Ql and W2. The estimates 
are highly correlated (p = 0.9). This high level of correlation is common to all pairs of channels and is caused by primary CMB fluctuations. 



